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, surface  wind  stress,  bottom  frictional  stress,  and  the  "inverted  barometer" 
effect,  but  omits  the  nonlinear  convective  acceleration  terms-. 

Surface  wind  stress  is  the  dominant  forcing  function  for  the  present 
application  of  the  numerical  model.  The  hurricane  wind  field  model  adopted  in 
this  study  conforms  to  the  1972  criteria  for  the  Standard  Project  Hurricane. 

An  ensemble  of  270  hypothetical  hurricanes  was  produced  by  compounding  three 
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forward  translation  speeds,  and  five  central  pressure  values.  The  central  pres¬ 
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Statistical  analyses  were  performed  for  three  beachfront  sites  in  the 
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open-coast  surge  and  for  the  net  setup  dug  to  the  peak  surge  and  wave  setup  due 
to  breaking.  For  each  of  the  nine  combiriatT&hs  of  landfall  site  and  beachfront 
site,  a  set  of  nomograms  was  produced.  These  provide  an  approximate  representa¬ 
tion  of  the  relation  between  peak  setup,  forward  speed,  radius  to  maximum  wind, 
and  central  pressure  deficit  for  each  of  the  three  track  headings.  The  nomo¬ 
grams  have  an  overall  root-mean-square  percentage  error  of  7.0  percent  when 
compared  with  individually  computed  peak  surge/wave  setup  values.  Among  the 
27  individual  nomograms,  the  percentage  error  covered  a  range  from  3.5  percent 
to  10. U  percent.  Maximum  absolute  error  in  the  worst  case  was  less  than  1  ft. 
This  level  of  accuracy  is  consistent  with  the  approximations  and  uncertainties 
inherent  in  the  numerical  models  used  in  the  study. 
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PREFACE 

A  request  for  a  computer  numerical  model  study  of  hurricane  surge 
setup  stage-frequency  relations  for  Dade  County,  Florida,  was  initiated 
by  the  District  Engineer,  U.  S.  Army  Engineer  District,  Jacksonville 
(SAJ),  in  March  1979.  Funds  for  the  U.  S.  Army  Engineer  Waterways  Ex¬ 
periment  Station  (WES)  to  conduct  the  study  were  authorized  via  DA  Form 
25UU,  Order  No.  08-123-ENG-0118-79,  dated  26  April  1979. 

The  model  study  was  conducted  during  the  period  June  1979  through 
January  1980  by  personnel  of  the  Wave  Processes  Branch  (WPB),  Wave 
Dynamics  Division  (WDD),  Hydraulics  Laboratory,  WES,  under  the  direction 
of  Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory;  Mr.  F.  A. 
Herrmann,  Jr.,  Assistant  Chief  of  the  Hydraulics  Laboratory;  Dr.  R.  W. 
Whalin,  Chief  of  the  WDD;  and  Mr.  C.  E.  Chatham,  Jr. ,  Chief  of  the  WPB. 
The  investigation  was  supervised  and  executed  by  Dr.  C.  E.  Abel,  Research 
Oceanographer,  WDD.  This  report  was  prepared  by  Dr.  Abel. 

During  the  course  of  the  investigation,  liaison  between  SAJ  and 
WES  was  maintained  by  means  of  telephone  communications  and  periodic 
progress  reports.  Technical  liaison  with  SAJ  was  provided  through  the 
services  of  Mr.  A.  0.  Hobbs. 

Commanders  and  Directors  of  WES  during  the  conduct  of  this  inves¬ 
tigation  and  the  preparation  and  publication  of  this  report  were 
COL  John  L.  Cannon,  CE,  and  COL  Nelson  P.  Conover,  CE.  Technical  Di¬ 
rector  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 
verted  to  metric  (Si)  units  as  follows: 


_ Multiply _ 

feet 

inches  of  Hg  (32°  F) 
knots  (international) 
miles  (U.  S.  nautical) 
miles  (U.  S.  statute) 


By _ 

0.301+8 

3386.38 

0.511*1+1+1+1+1+ 
1.852 
1. 60931+1+ 


To  Obtain 
metres 
pascals 

metres  per  second 

kilometres 

kilometres 
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HURRICANE  SURGE  STAGE-FREQUENCY  ANALYSIS 
FOR  DADE  COUNTY,  FLORIDA 


PART  I :  INTRODUCTION 


1.  Since  the  occurrence  of  the  disastrous  hurricane  of  September 
1926,  at  least  12  hurricanes  and  U  tropical  storms  have  made  landfall 
within  50  to  60  nautical  miles*  of  the  greater  metropolitan  Miami  area 
(Neumann  et  al.  1978).  In  addition  to  these  storms  which  entered  land 
from  the  Atlantic  Ocean,  a  number  of  other  storms  whose  tracks  passed 
beyond  the  60-mile  radius  or  which  crossed  the  Florida  peninsula  and 
subsequently  exited  into  the  Atlantic  have  produced  high  storm  tides  in 
the  Miami  area. 

2.  In  response  to  the  enormous  property  damage  caused  by  the  1926 
storm,  a  continuing  effort  has  been  made  to  minimize  the  threat  of  storm 
surge  damage  through  construction  of  beach  erosion  and  hurricane  protec¬ 
tion  structures.  The  most  recent  series  of  improvements  are  being  under¬ 
taken  by  the  U.  S.  Army  Engineer  District,  Jacksonville  (SAJ).  The  im¬ 
provements  apply  to  a  10. 5-mile  segment  of  the  ocean  coast  of  Dade 
County,  extending  from  Government  Cut  to  the  northern  end  of  Haulover 
Beach  Park.  The  project  provides  for  beach  erosion  control  and  hurricane 
protection  along  the  9.3  miles  of  shoreline  extending  northward  from 
Government  Cut  to  Bakers  Haulover  Inlet  by  the  creation  of  a  protective 
hurricane  surge  dune  with  a  20- ft  crown  at  11.5  ft  above  mean  low  water 
(mlw)  and  side  slopes  of  1  on  5,  down  to  a  protective  beach,  with  a 
level  berm  50  ft  wide  at  el  9  ft  above  mlw,  and  a  natural  slope  seaward. 
For  the  1.2  miles  of  Haulover  Beach  Park,  beach  erosion  control  measures 
are  being  provided  by  the  creation  of  a  protective  beach  of  the  same 
dimensions  as  previously  described.  In  addition,  studies  are  currently 
being  conducted  by  SAJ  to  determine  the  feasibility  of  extending  these 


*  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  presented  on  page  3. 
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improvements  northward  to  the  north  county  line. 

3.  The  numerical ■ study  described  in  this  report  was  undertaken  to 
examine  the  probability  of  occurrence  of  hurricane-induced  water-level 
anomalies.  For  purposes  of  this  study,  an  idealized  beach  representing 
the  beach  erosion  control  and  hurricane  protection  features  of  the  proj¬ 
ect  improvements  was  assumed  to  extend  from  Government  Cut  to  the  north¬ 
ern  boundary  of  Dade  County.  Figure  1  shows  the  area  of  interest  and 
also  indicates  locations  at  which  water  level  is  computed  by  the  numeri¬ 
cal  models  employed  in  this  study.  Two  main  kinds  of  data  were  obtained 
from  this  study:  (a)  stage -frequency  curves  for  peak  surge  and  combined 
storm  surge-wave  setup  peak  values  at  three  beachfront  sites  and  (b)  nomo¬ 
grams  providing  for  approximate  calculation  of  surge  and  setup  for  given 
hurricane  parameters  and  landfall  site.  An  ensemble  of  270  hypothetical 
hurricanes  was  used  to  deduce  the  probability  distributions  of  combined 
peak  surge  and  wave  setup,  and  to  construct  the  prediction  nomograms. 

1*.  Results  of  this  study  may  be  used  to  estimate  the  probability 
of  occurrence  of  storm  tides  in  the  area  of  interest,  assuming  the  beach 
has  been  uniformly  modified  to  match  the  idealized  geometry  described 
later.  The  validity  of  the  results  is  also  subject  to  the  limitations 
on  the  ensemble  of  hypothetical  hurricanes.  No  extreme  events  were  in¬ 
cluded  in  this  study,  beyond  the  100-year  storm  (based  on  frequency  of 
occurrence  of  central  pressure  for  southeast  Florida).  Similarly,  re¬ 
sults  of  this  study  provide  a  scheme  for  estimation  of  peak  surge  and 
wave  setup  for  storms  contained  in  the  given  ensemble.  Extrapolation 
to  storms  not  in  the  study  ensemble  is  not  recommended.  However,  inter¬ 
polation  over  parameters  of  the  ensemble  storms  may  be  attempted,  but 
the  results  should  be  considered  approximate  since  the  intrinsic  pro¬ 
cesses  involved  are  not  necessarily  linear. 
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PART  II:  COMPUTATIONAL  TECHNIQUES 


Surge  Computations 

5.  The  basic  numerical  model  used  for  surge  computations  in  this 
study  is  the  two-dimensional  storm  surge  computer  code,  SSURGE  III.  The 
model  solves  an  appropriate  system  of  three  coupled  partial  differential 
equations.  Two  equations  govern  the  evolution  of  the  components  of 
vertically  integrated  volume  transport,  and  the  third  equation  describes 
the  condition  for  continuity  of  the  transport.  A  complete  description 
of  development  of  the  model  and  a  survey  of  problems  to  which  it  has 
previously  been  applied  are  given  in  Wanstrath  et  al.  (1976)  and 
Wanstrath  (1977).  A  concise  summary  of  the  major  features  of  this  model 
follows . 

6.  The  more  traditional  approach  of  applying  the  equations  of 
motion  in  discrete  form  to  a  rectangular  coordinate  system  is  not  fol¬ 
lowed  by  the  SSURGE  III  code.  Recognizing  that  most  of  the  Atlantic 
and  Gulf  coasts  of  the  United  States  are  better  described  by  smoothly 
varying  curves  than  by  the  stairstep  outline  inherent  in  a  rectangular 
grid,  this  model  uses  the  shoreline  itself  to  define  the  appropriate 
coordinate  system.  This  choice  of  natural  coordinates  leads  to  formu¬ 
lation  of  the  model  in  an  orthogonal  curvilinear  coordinate  system.  A 
region  bounded  by  natural  curved  boundaries  is  mapped  into  a  rectangular 
region  of  a  uniquely  defined  image  space.  Grid  line  spacings  in  this 
rectangular  region  are  then  subjected  to  an  additional  mapping  to  allow 
for  high  spatial  resolution  in  regions  of  interest  and  lower  resolution 
in  regions  of  the  grid  at  which  detailed  information  is  not  required. 

The  computational  grid  for  this  study  is  shown  in  Figure  2,  which  gives 
locations  at  which  water-surface  elevation  is  computed.  The  two  com¬ 
ponents  of  integrated  transport  are  computed  at  midpoints  of  the  edges 
of  the  quadrilaterals  defined  by  these  elevation  points. 

7*  In  the  doubly  transformed  coordinate  system,  the  continuous 


form  of  the  equations*  of  motion  is  given  by  the  following. 


3t  fQT  +  yF  3S  TS  “  °S  ^ 

Tt-  +  +  vF  3T  ^  “  tT  "  °T  ^ 

f  f^[uls  <F«S>  <fv]*°  (3> 


where 

Q  =  volume  transport  per  unit  width 
S  =  shoreline  coordinate 
t  =  time 

f  =  Coriolis  parameter 
T  =  coordinate  normal  to  S 
g  =  acceleration  due  to  gravity 
D  =  total  water  depth 

F,y,v  =  variable  scale  factors  describing  the  combined  effects  of 
the  conformal  and  stretching  transformations  applied  to 
the  original  curvilinear  coordinate  system 

H  =  sea  surface  elevation  relative  to  a  specified  datum 

Hg  =  hydrostatic  sea  surface  elevation  due  to  atmospheric 
pressure  anomaly  ("inverted  barometer"  effect) 

x  =  wind  stress  divided  by  water  density 

a  =  bottom  stress  (friction)  divided  by  water  density 

Equations  1  through  3  are  then  discretized  and  the  resulting  system  of 

coupled  algebraic  equations  are  solved  by  an  explicit  time-stepping 

algorithm. 

8.  Boundary  conditions  for  the  computational  region  are  (a)  in¬ 
verted  barometer  height  specified  at  the  seaward  boundary,  (b)  "leaky 
wall"  condition  at  the  shore  (concept  described  in  Wanstrath  1978), 

*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  in  the 
Notation  (Appendix  A). 


9 


and  (c)  vanishing  normal  gradient  of  normal  component  of  transport 
(i.e.,  3Qg/3S  =  0)  at  the  two  lateral  faces.  For  this  particular  study, 
the  usual  interior  initial  condition  of  no  motion  was  replaced  by  an 
initial  condition  conducive  to  smoothing  the  solution  for  the  early 
time-steps.  An  approximate  balance  between  Coriolis  acceleration  and 
the  effects  of  horizontal  pressure  gradient  was  assumed.  That  is. 


fo  (U) 

pF  3S  K  1 

f QS  =  vF  3T~  ^ 5  ^ 


Equations  U  and  5  ignore  the  local  acceleration  and  the  wind  and  bottom 
frictional  stress  terms.  The  SSURGE  III  code  also  takes  no  account  of 
permanent  current  systems  such  as  the  Florida  Current  segment  of  the 
Gulf  Stream  which  flows  offshore  of  the  study  region.  Wave-current 
interactions  are  essentially  nonlinear  phenomena  and  cannot  be  adequately 
treated  by  the  simplified  equations  of  motion  utilized  in  this  model. 
Approximate  values  of  initial  transport  are  thus  provided  by  simply 
solving  these  equations,  given  the  initial  distribution  of  the  inverted 
barometer  height  HD  .  To  complete  the  specification  of  initial  condi- 

.D 

tions,  the  surface  elevation  was  equated  to  the  inverted  barometer 
height.  These  two  conditions  are  not  entirely  accurate,  but  they  were 
found  effective  in  eliminating  transients  from  the  solution  for  early 
time-steps. 

9.  The  final  characteristic  of  the  computational  grid  to  be  dis¬ 
cussed  is  specification  of  bottom  topography.  The  shelf  margin  is  un¬ 
usually  close  to  the  coast  in  the  area  of  interest.  In  cases  of  broad, 
slowly  varying  shelf  margins,  the  seaward  boundary  is  typically  taken 
out  to  the  300-  or  600-ft  contour.  Offshore  from  Biscayne  Bay,  these 
depths  are  found  only  a  few  miles  seaward.  The  constriction  in  the  cur¬ 
vilinear  region  caused  by  adhering  to  the  600-ft  contour  was  found  to 
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produce  serious  instabilities  in  the  numerical  solution  of  the  equations 
of  motion.  Stability  was  achieved  by  using  an  offshore  boundary  ob¬ 
tained  by  moving  the  shoreline  boundary  parallel  to  itself  in  the  sea¬ 
ward  direction  a  distance  of  51  nautical  miles.  Moving  from  the  shore¬ 
line  in  the  seaward  direction,  actual  values  of  depth  were  used  until 
the  1200-ft  contour  was  reached  (NOS  Chart  llU60  was  used  to  obtain 
these  values).  Beyond  this  line,  a  constant  depth  of  1200  ft  was  ex¬ 
tended  to  the  seaward  boundary  of  the  grid.  Although  this  process  com¬ 
pletely  ignores  the  presence  of  the  Bahamas  Bank,  this  idealized  topog¬ 
raphy  was  essential  for  obtaining  stable  solutions  for  all  hypothetical 
storms  used  in  this  study. 

Wave  Setup  Computations 

10.  In  addition  to  determination  of  the  peak  surge  at  the  coast, 
the  increment  of  water  level  due  to  the  setup  of  breaking  storm  waves 
was  considered.  The  amplitude  H  and  the  period  T  of  the  signifi- 

D  D 

cant  storm  wave  were  determined  from  hurricane  model  parameters,  follow¬ 
ing  the  procedure  detailed  in  paragraph  27.  This  significant  wave  was 
carried  landward  along  approximate  orthogonals  until  breaking  occurred. 
Characteristics  of  the  breaking  wave  were  used  to  estimate  wave  setup 
on  the  idealized  beach  representing  the  projected  improvements  in  the 
area  of  interest.  Details  of  this  procedure  follow. 

11.  To  assure  that  true  deepwater  conditions  existed  for  the  sig¬ 
nificant  storm  wave,  the  seaward  end  of  the  orthogonals  was  taken  at  the 
800-ft  contour  (located  by  hand-contouring  the  depth  values  on  NOS 
Chart  1117*0 •  Each  orthogonal  was  divided  into  four  segments,  from  the 
800-ft  contour,  shoreward  to  the  mlw  line.  The  scheme  used  to  select 
these  segments  was:  (a)  800-ft  contour,  shoreward  to  a  point  12,000  ft 
offshore;  (b)  12,000  ft  offshore  to  a  variable  point  selected  from  NOS 
Chart  SC-11467  to  idealize  nearshore  topography;  (c)  variable  point  to 
15-ft  contour  (assumed  to  be  600  ft  offshore  from  mlw  line);  and 

(d)  15-ft  contour  to  mlw  line.  Landward  of  the  mlw  line  an  arbitrarily 
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wide  beach  slope  of  1:20  was  assumed.  This  sloping  beach  serves  as  a 
backdrop  against  which  the  combined  effects  of  storm  surge  and  wave 
setup  are  projected.  Total  stage  values  may  accordingly  be  interpreted 
as  "overtopping"  if  a  specific  elevation  above  mlw  is  used  as  a  cutoff 
valve. 

12.  The  idealized  topography  described  above  exists  in  the  ab¬ 
sence  of  surge  setup.  If  a  surge  S  was  computed  at  the  coast,  it  was 
assigned  a  location  at  the  mlw  line.  Locations  on  the  orthogonals  were 
assigned  augmented  total  depths  by  adding  zero  at  the  800-ft  contour  and 
using  linear  interpolation  for  points  between  mlw  and  the  offshore 
limit.  These  augmented  depths  will  be  referred  to  by  the  symbol  d  . 
Since  there  is  now  a  positive  water  elevation  at  the  mlw  line,  there 
must  be  water  on  the  sloping  beach  beyond  this  line.  This  is  provided 
by  projecting  the  peak  surge  S  until  the  sloping  beach  is  encountered. 
For  the  wave  setup  study,  there  is  now  a  fifth  segment  to  be  considered. 
Orthogonals  were  defined  for  the  nine  gage  locations  shown  in  Figure  1. 
Augmented  depths  were  used  to  solve  the  wave  modification  problem  for 
the  significant  storm  wave. 

13.  Let  Hg  and  Tg  be  the  height  and  period,  respectively,  of 
the  significant  storm-induced  wave.  These  values  are  taken  to  exist  at 
the  800-ft  contour.  Wave  modification  is  assumed  to  occur  along  orthog¬ 
onals,  so  refraction  effects  are  neglected.  Bottom  percolation  and 
nonlinear  wave-wave  interactions  also  are  neglected.  Thus  only  the 
effects  of  geometric  spreading  (shoaling)  and  bottom  friction  are  con¬ 
sidered.  The  method  of  computation  followed  is  due  to  Bretschneider  and 
Reid  (195M. 

lU.  Consider  the  ith  segment  of  an  orthogonal.  Its  length  is 
Ax^  and  the  water  depths  are  d^  on  the  seaward  end  and  d^+1  on  the 
landward  end.  If  is  the  wave  height  on  the  seaward  end  of  the 

segment  and  H^+1  is  the  modified  wave  height  on  the  landward  end,  then 

Hi+1  =  Ky  H.  (6) 
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The  terms  multiplying  Ax^  in  Equation  7  are  defined  as  follows.  Th 
factor  A  is  dependent  on  bottom  friction. 


where  A  is  a  friction  factor,  taken  as  0.012  in  this  study.  The 
factor  Kc  is  the  linear  long-wave  form  of  the  shoaling  factor,  com- 

O 

puted  here  as 


KS  = 


^  sinh  Urr(d*/L?) 

tanh  2ir(d*/L*)  Uir(d?/L*)  +  sinh  U7i(d 


•)  l1/2 

r _ 

ir(d*/L*5" 


where  d?  =  (di  +  di+^)/2  and  the  wavelength  L*  is  the  solution  to 
the  linear  long-wave  characteristic  relation 


Finally,  the  third  factor,  a  geometric  form  factor,  is  given  by 


* .  ^  r  ks  i3 

*F  ^g2  I  sinh  2ir(d*/L*)J 

The  shoaling  wave  is  assumed  to  have  broken  in  the  ith  segment  if 


( 


the  condition  _>  0.78  d^+1  is  met.  The  shoaling  computations  are 

started  by  setting  H  =  Hc  . 

■L  o  t  h 

15.  If  the  shoaling  wave  has  broken  in  the  i  segment,  then 

linear  interpolation  is  used  to  determine  the  location  and  the 

depth  d^  of  breaking.  Following  a  recommendation  from  the  Shore 

Protection  Manual  (CERC  1977 ) »  the  net  wave  setup  on  a  normal  sloping 

beach  is  taken  to  be  represented  by  two  components 


After  computation  of  has  been  completed,  the  total  storm-induced 

stage,  due  to  peak  surge  and  wave  setup,  is  given  by  HgW  =  S  +  . 

The  quantity  Hsw  is  the  primary  variable  considered  in  the  subsequent 
stage-frequency  analysis.  Interpolation  between  H  and  the  wave 

Dn 

height  in  the  breaking  segment  provides  an  approximate  means  of  examin¬ 
ing  water  elevation  in  the  nearshore  zone.  In  particular,  an  interpola¬ 
tion  procedure  of  this  type  was  used  to  determine  wave-height  distribu¬ 
tions  at  the  position  of  the  still-water  15-ft  contour.  Results  of  this 
computation  are  presented  in  Plates  1-3.  Each  of  these  distributions 
is  based  on  a  sample  of  810  estimated  heights,  comprising  270  values  at 
each  of  the  three  model  gages  in  the  beachfront  segment. 
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PART  III:  HURRICANE  WIND  FIELD  MODEL 


Wind  Field  Model  Selection 


16.  Surface  wind  stress  is  the  dominant  driving  force  in  generat¬ 
ing  the  peak  surge  impacting  the  study  area.  Accordingly,  a  recognized 
and  accepted  model  for  the  hurricane  winds  should  be  used.  This  crite¬ 
rion  is  satisfied  by  the  Standard  Project  Hurricane  (SPH).  The  revised 
criteria  published  by  the  National  Weather  Service  in  1972  form  the 
basis  for  computational  aspects  of  the  wind  field  model  (NOAA  1972). 
Implementation  of  these  criteria  was  accomplished  with  a  computer  code 
designed  for  use  in  surge  studies  (Reid  1979)* 

17.  The  SPH  wind  model,  as  defined  by  certain  geometric  relations, 
describes  the  distribution  of  wind  velocity  for  overwater  conditions. 

The  magnitude  of  the  wind  depends  mainly  on  radius  from  the  low  pressure 
center  to  the  point  of  maximum  wind  speed  (R^),  the  magnitude  of  the 
central  pressure  deficit  relative  to  an  assumed  constant  peripheral 
pressure  (dP),  and  forward  translational  speed  of  the  storm  (VT).  The 
computer  code  used  in  this  study  adheres  to  the  1972  criteria  for 
specification  of  an  overwater  wind  field.  Modification  of  the  wind 
field  by  the  presence  of  land  also  is  included  in  the  model  through 
introduction  of  fetch-dependent  reduction  factors  for  wind  speed.  The 
procedure  recommended  by  the  1972  criteria  for  computation  of  the  re¬ 
duction  factors  is  included  in  the  SPH  computer  code.  Fetch  dependence 
of  wind  speed  near  the  coastal  zone  thus  adds  storm  track  heading  (0T) 

to  the  list  of  pertinent  factors  for  characterization  of  the  wind  field. 
This  single  parameter  is  sufficient  for  specification  of  the  track 
shape  since  the  SPH  model  requires  linear  tracks  for  storm  motion. 

18.  Ease  of  application  and  small  number  of  input  parameters  were 
significant  in  selecting  the  SPH  wind  field  model  since  a  fairly  large 
number  (270)  of  storms  w«“r**  used  in  the  surge  stage-frequency  analysis. 
The  SPH  model,  with  inclusion  of  fetch-dependent  wind-reduction  factors, 
appears  more  than  adequate  for  purposes  of  a  stage-frequency  analysis. 

It  is  easy  to  implement,  relatively  inexpensive  to  operate,  and  gives 
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a  reasonable  representation  of  hurricane  wind  fields,  including  distor¬ 
tion  of  the  wind  field  as  the  storm  approaches  shore. 

Selection  of  Storm  Parameters 


19.  The  first  set  of  parameters  to  be  considered  are  those  which 
specify  the  storm  track.  Since  SPH  tracks  are  linear,  only  a  single 
angle  is  required  for  orientation  of  the  track.  Absolute  location  of 
the  track  is  fixed  by  specifying  the  landfall  site  on  the  coast.  Three 
landfall  sites  were  selected  for  this  study.  They  are  marked  in  Fig¬ 
ure  2  with  circled  dots  and  the  labels  LF^,  LF^y  and  LF^.  The  first 
location  selected  was  LF^,  near  the  area  of  Atlantic  Heights.  The 
rationale  is  the  known  hurricane  characteristic  of  maximum  winds  lying 
in  the  quadrant  to  the  right  of  the  advancing  storm.  Landfall  site 
LF^  can  reasonably  be  expected  to  be  critical  for  maximum  surge  and 
wave  setup  in  the  area  of  interest.  Landfall  sites  LF^  and  LF^  were 
selected  an  equal  distance  upcoast  and  downcoast  from  the  critical  land¬ 
fall  site,  but  still  near  enough  the  area  of  interest  to  produce  signifl 
cant  effects.  The  possible  range  of  storm  track  headings  was  examined 
in  the  series  of  hurricane  tracks  published  by  Neumann  et  al.  (1978). 
These  tracks  revealed  that  the  majority  of  storms  making  landfall  within 
60  miles  of  the  area  of  interest  were  confined  to  headings  in  the  sector 
0,p  ■  120  deg  true  to  0,^  =  170  deg  true.  Track  heading  is  used  here 
in  the  meteorological  sense;  i.e.,  true  compass  direction  from  which 

the  storm  is  approaching  the  landfall  site.  Three  track  headings  were 
selected  for  this  study:  120  deg,  ll*5  deg,  and  170  deg  true.  To  com¬ 
plete  specification  of  the  tracks,  all  hypothetical  storms  were  modeled 
for  a  total  distance  of  160  nautical  miles,  from  start  of  simulation  to 
time  of  landfall.  Total  time  for  each  simulation  of  storm  movement  was 
fixed  at  2h  hr,  thus  assuring  that  all  hypothetical  storms  were  well 
past  the  shoreline  at  the  end  of  the  simulation. 

20.  The  other  three  important  parameters  (R^,  V,p,  AP)  were  se¬ 
lected  by  consideration  of  historical/statistical  data.  Both  radius 
to  maximum  wind  and  storm  translational  speed  were  obtained  from  the 
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original  data  base  used  in  defining  the  SPH  (Graham  and  Nunn  1959). 

This  publication  lists  ranges  of  R^  and  V^,  values  appropriate  to  the 
southeastern  coast  of  Florida,  including  the  area  of  interest.  Dis¬ 
counting  extreme  values,  two  values  of  R^  and  three  values  of  VT 
were  selected.  The  radius  to  maximum  wind  was  chosen  to  have  the  values 
8  and  15  nautical  miles.  Similarly  discarding  extreme  values,  the 
values  for  translational  speed  selected  were  8,  12,  and  16  knots. 

21.  The  statistical  distribution  of  central  pressure,  in  inches 
of  Hg,  for  the  southeastern  Florida  coast  is  provided  in  a  report  by  SAJ 
(Jacksonville  District  1963).  Data  used  in  deriving  this  distribution 
cover  hurricanes  occurring  between  the  years  1900  and  i960.  The  values 
of  central  pressure,  in  inches  of  Hg,  selected  for  inclusion  in  this 
study  were  28.33,  27-7^.  27. Ul,  26.95>  and  26.78.  Using  the  indicated 
pressure-frequency  relation,  it  is  easily  shown  that  these  values  of 
central  pressure  correspond  to  mean  recurrence  intervals  of  5,  10,  20, 

50,  and  100  yr,  respectively.  To  obtain  values  of  AP  ,  the  peripheral 
pressure  was  assumed  to  be  one  standard  atmosphere,  29.92  inches  of  Hg. 

22.  It  has  probably  been  noticed  that  no  mention  has  yet  been 

made  of  maximum  wind  speed  as  an  important  parameter  in  determining 

surge  storm-wave  characteristics.  The  maximum  wind  speed  (V.,)  is  indeed 

M 

important  but  is  not  a  free  parameter  in  the  SPH  model;  rather,  it  is  a 
function  of  three  parameters.  The  form  of  this  dependence  is 

VM  =  6U.6(AP)1/2  -  0.56  Rj^  +  0.5  vT  (15) 

where  a  latitude  of  26  deg  N  has  been  assumed  in  computing  the  coeffi¬ 
cients  for  this  equation.  The  same  latitude  was  used  to  determine  the 
constant  Coriolis  parameter  required  by  the  numerical  surge  model. 

23.  Recapitulating  the  selection  of  parameters,  the  ensemble  is 
described  by:  three  landfall  sites,  three  track  headings,  two  radii  to 
maximum  wind,  three  translational  speeds,  and  five  central  pressure  defi¬ 
cits.  Thus  the  ensemble  of  different  hypothetical  hurricanes  used  in 
this  study  has  270  members.  Organization  and  analysis  of  computational 
results  for  these  storms  are  discussed  in  the  next  part  of  this  report. 
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PART  IV:  STAGE-FREQUENCY  ANALYSIS 

Peak  Surge  Statistics 

2k.  The  SSURGE  III  program  was  run  for  each  member  of  the  270 
hypothetical  hurricane  ensemble.  Maximum  water  elevation  was  recorded 
at  each  of  the  nine  numerical  gages  in  the  area  of  interest.  In  order 
to  increase  the  sample  size  available  for  statistical  analysis  at  a  site, 
the  gages  were  grouped  by  threes  to  represent  observations  at  the  beach¬ 
front  sites  labeled  BF^,  BF2,  and  BF^  in  Figure  1.  The  data  comprising 
peak  surge  values  contain  samples  from  three  distinct  populations,  deter¬ 
mined  by  location  of  the  landfall  site.  Accordingly,  landfall  site  is 
used  as  a  conditional  parameter  when  any  statistical  assertions  are  made 
about  the  sample  values.  That  is,  any  assertion  should  have  a  qualifying 
phrase  such  as,  "given  landfall  site  X  ,  then  the  probability  of  the 
peak  surge  exceeding  Y  feet  at  beach  front  Z  is  75  percent." 

25.  In  light  of  this  use  of  landfall  site  to  separate  the  peak 

surge  values  into  three  populations,  it  should  be  apparent  that  each 
beachfront  site  is  represented  by  a  sample  of  270  peak  surge  values.  Es¬ 
timates  of  the  cumulative  distribution  function  (P)  for  peak  surge  at 
each  beachfront  site  were  computed  with  respect  to  a  sampling  window 
width  of  0.5  ft.  A  set  of  values  =  0. 5  n  is  generated  for  n  =  1, 

2,  ...»  N  ,  where  N  is  the  smallest  integer  such  that  exceeds 

the  largest  peak  surge  value  in  a  given  sample.  Then  the  N  numbers, 

F  ,  are  determined  such  that  F  is  the  number  of  sample  values  less 
n  n 

than  or  equal  to  .  Estimates  of  P  are  then  given  by 

F 

P(W  )  =  Prob(w  <  W  )  =  rr2-  (16) 

n  —  n  M 

where  M  =  270  is  the  sample  size  and  w  is  the  intrinsic  random 
variable  (peak  surge)  from  which  the  computed  values  are  taken  to  be 
a  set  of  independent  random  samples. 

26.  The  customary  form  for  presentation  of  the  cumulative 
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distribution  function  for  peak  surge  is  in  terms  of  its  complement,  the 
exceedance  (Q),  defined  by 

Q(W  )  =  Prob(w  >  w  )  =  1  -  P(W  )  (17) 

That  is,  the  fundamental  information  of  interest  in  discussing  peak 
surge  distributions  is  the  probability  of  exceeding  a  specified  stage 
level.  Accordingly,  curves  of  Q  are  presented  in  Plates  k-6.  Each 
plate  comprises  three  exceedance  curves,  each  referring  to  a  landfall 
site,  as  discussed  in  the  previous  paragraph.  An  alternative  method  of 
presenting  the  peak  surge  distribution  is  to  plot  peak  surge  versus 
recurrence  interval  of  the  central  pressure  deficit.  In  the  present 
case,  there  are  five  values  for  central  pressure  deficit,  with  recurrence 
intervals  as  described  in  paragraph  21.  For  each  beachfront  site  there 
is  a  sample  of  162  peak  surge  values  for  each  of  the  five  recurrence 
intervals.  These  samples  were  used  to  compute  the  curves  given  in 
Plates  7-9-  Each  plate  shows  three  curves — a  mean  peak  surge  curve 
(with  error  bounds  indicated  at  computational  points),  and  minimum/ 
maximum  peak  surge  curves  to  illustrate  the  range  of  values  encountered 
at  the  beachfront  site. 


Combined Peak  Surge/Wave  Setup  Statistics 


27 .  The  peak  surge  values  recorded  at  the  beachfront  gages  were 
used  to  compute  wave  setup  associated  with  the  significant  waves  for 
each  hypothetical  hurricane.  As  described  in  PART  II,  the  computational 

scheme  requires  for  input  the  height  (H  )  and  the  period  (Te)  of  the 

b  b 

deepwater  significant  storm  wave.  Both  these  wave  characteristics  may 
be  estimated  in  terms  of  basic  hurricane  parameters  following  a  method 
outlined  in  the  Shore  Protection  Manual  (CERC  1977).  The  pertinent 
formulas  are 
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(19) 


Tg  =  8.6  exp 


\  200  / 


1  + 


rO.lOU  V„ 


^ M 


where  ^  »  AP  ,  are,  respectively,  radius  to  maximum  wind  (nautical 

miles),  central  pressure  deficit  (inches  of  Hg),  and  forward  speed 

(knots).  The  maximum  wind  speed.  (V,.  ,  knots)  is  given  by  the  usual  SPH 

M 

formula  (Equation  15).  Within  the  storm  ensemble  there  are  30  values 
for  both  H0  and  T  .  Frequency  distributions  of  these  deepwater 

D  O 

storm-wave  characteristics  are  illustrated  in  Plate  10.  The  small  sample 
size  necessitated  use  of  fairly  broad  sampling  intervals — 0.5  sec  for 
the  Tg  distribution  and  2.0  ft  for  the  distribution.  Since  linear 

long-wave  theory  is  used  for  the  wave  setup  computations,  the  frequency 
distribution  of  wave  period  in  the  nearshore  region  may  be  assumed  to  be 
identical  to  the  deepwater  Tg  distribution. 

28.  Calculations  were  performed  at  each  gage,  for  each  of  the 
270  hypothetical  storms  to  produce  the  required  2^30  estimates  of  Hgw  , 
the  combined  peak  surge  and  wave  setup  at  the  beach.  As  before,  results 
were  separated  into  the  three  basic  populations,  and  then  grouped  into 
three  beachfront  samples  of  270  values.  Exceedance  probabilities  were 
computed  and  results  from  this  computation  are  presented  in  Plates  11-13. 
An  alternative  presentation  of  these  peak  setup  results  in  terms  of 
central  pressure  deficit  recurrence  interval  is  shown  in  Plates  lk-l6. 

The  format  is  identical  to  that  of  the  peak  surge-recurrence  interval 
curves  described  in  paragraph  26. 

Prediction  Nomograms  for  Peak  Setup 

29.  The  two  types  of  exceedance  distributions  described  above  are 
characterizations  of  storm  populations  whose  individual  members  are 
functions  of  a  number  of  assumed  independent  parameters —  9^  , 

AP  ,  and  .  Details  of  individual  storms  are  more  or  less  lost  in 
calculating  population  statistics.  It  is  possible,  however,  to  use  the 
individual  storm  parameters  to  approach  the  problem  of  predicting  peak 
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setup  at  each  of  the  beachfront  sites.  As  used  here,  "prediction"  means 
obtaining  a  good  approximation  to  calculated  Hgw  values  given  the 
above  four  parameters  as  input  data.  Such  an  approximate  prediction 
technique  is  derived  for  the  results  of  this  study. 

30.  Again,  landfall  site  is  used  as  a  conditional  parameter.  Of 

the  four  input  parameters,  0^  has  no  direct  dynamical  relation  to 
Hgy  .  It  does  have  a  qualitative  impact  on  the  results  since  it  speci¬ 
fies  the  "crossing  class"  of  a  storm.  Accordingly,  0^,  is  also  used 
as  a  conditional  parameter,  leaving  ,  AP  ,  and  to  be  directly 

related  to  peak  setup.  For  a  specified  landfall  site,  and  given  a  par¬ 
ticular  storm  track  heading,  the  beachfront  value  of  Hgw  for  any  com¬ 
bination  of  other  parameters  was  taken  as  the  average  of  peak  setup  at 
the  three  gage  sites  defining  the  beachfront  site.  Repeated  examination 
of  the  behavior  of  Hgw  with  respect  to  ,  AP  ,  and  VT  revealed 
that  the  curve  of  Hg^  /V»T  versus  represented  all  computed 

results  with  an  overall  root-mean- square  (RMS)  percentage  error  of  7-0 
percent.  For  any  particular  landfall  site  and  beachfront  site,  the 
value  of  the  range  of  individual  RMS  percentage  error  ranged  from  3.5  to 
10. U  percent.  The  nomograms  representing  these  collective  results  are 
given  in  Plates  17-25.  There  is  a  family  of  nomograms  for  each  of  the 
nine  combinations  of  landfall  site  and  beachfront  site. 

31.  Given  the  levels  of  uncertainty  and  inaccuracy  involved  in 

the  computation  of  Hg^  ,  the  indicated  level  of  error  is  considered 
acceptable.  The  RMS  percentage  error  applicable  to  the  27  individual 
prediction  nomograms  is  shown  in  the  tabulation  below.  It  should  be 
noted  that  the  worst  case  error  in  fitting  Hgw/  V*T  versus  in¬ 

volves  an  RMS  estimated  error  of  0.9  ft  for  a  predicted  H  of  7.6  ft. 

bW 

The  entries  in  the  tabulation  may  be  used  to  provide  a  refined  estimate 
of  error,  or  the  blanket  figure  of  7-0  percent  may  be  adopted  for  the 
entire  ensemble  of  nomograms. 


Percentage  Error  for  Individual  Nomograms 
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32.  Representation  of  calculated  peak  setup  values  in  the  form  of 
nomograms  was  motivated  by  needs  of  convenience,  and  not  on  theoretical 
grounds.  The  alternative  to  the  nomogram  representation  is  to  con¬ 
struct  a  five-dimensional  table  for  each  of  the  three  selected  landfall 
sites.  The  particular  characterization  selected  was  used  because  it 
provided  a  sufficiently  accurate  means  of  estimating  peak  setup.  It 
should  be  noted  that  this  particular  parameterization  is  not  very  accu¬ 
rate  if  peak  surge  (S)  is  substituted  for  peak  setup  (Hg^).  This  possi¬ 
bility  was  examined,  with  selected  computations  showing  more  than  twice 
the  percentage  error  than  the  nomograms  included  here.  This  suggests 
that  the  contribution  to  Hgw  from  wave  breaking  is  dominant  in  deter¬ 
mining  the  overall  character  of  the  peak  setup  response  function.  Since 
the  deepwater  wave  parameters  depend  on  ,  AP  ,  and  ,  this  may 
account  for  the  success  of  the  nomogram  representation  selected.  There 
is,  of  course,  the  possibility  that  the  observed  accuracy  is  an  artifact 
of  the  geometry,  wind  field  model,  and  computational  techniques  used 

in  this  study.  No  inference  is  intended  with  respect  to  the  applicabil-  t 

ity  of  this  representation  to  stage-frequency  analyses  for  other  loca¬ 
tions,  or  utilizing  different  wind  field  or  hydrodynamic  models. 

33.  It  is  also  necessary  to  consider  the  limits  of  applicability  ^ 

of  computed  results  within  the  study  ensemble  of  hypothetical  hurricanes. 

The  ranges  of  all  pertinent  hurricane  parameters  were  selected  to  avoid 
extremal  events.  This  fact  is  particularly  evident  in  selection  of 
central  pressure  values.  The  computed  occurrence  probabilities  used  were 
based  on  a  sample  of  approximately  20  storms  occurring  over  a  60-yr 


22 


interval.  This  sample  size  is  too  small  for  accurate  determination  of 
extremal  probabilities.  Thus,  the  computed  peak  surge  and  peak  setup 
exceedance  curves  do  not  include  any  pressure  events  with  a  mean  recur¬ 
rence  interval  exceeding  100  yr.  This  fact  should  be  kept  in  mind  when 
one  is  examining  and  interpreting  the  estimated  exceedance  functions. 

3b.  A  similar  restriction  applies  to  storm  track  heading.  Only 

three  distinct  values  were  used.  There  is  no  obvious  trend  in  H0„ 

SW 

response  with  respect  to  0,^  variation,  indicating  that  interpolation 
or  extrapolation  to  other  crossing  angles  may  produce  uncertain  values. 
This  is  also  the  case  if  landfall  sites  other  than  the  three  selected 
sites  are  examined.  Extensions  to  parameters  outside  the  ranges  se¬ 
lected  for  this  study  may  yield  reasonable  results ,  but  no  estimate  of 
their  accuracy  can  be  provided.  The  preferred  method  of  using  the  com¬ 
puted  nomograms  is,  for  any  set  of  storm  parameters,  to  select  values 
closest  to  a  particular  value  included  in  this  study.  Accordingly,  for 
parameter  values  well  outside  the  ranges  included  in  this  study,  the 
nomograms  should  not  be  expected  to  yield  peak  setup  values  consistent 
with  an  independent  computation  with  the  desired  parameters. 


PART  V:  SUMMARY 


3S  A  methodology  for  hurricane  surge  stage-frequency  analysis 
for  Dade  County,  Florida,  has  been  presented.  Computational  aspects  of 
this  methodology  include  appropriate  numerical  representations  of  the 
basic  hydrodynamics  and  of  the  hurricane  wind  field.  The  vertically 
integrated  hydrodynamics  model  was  based  on  an  orthogonal  curvilinear 
coordinate  system,  with  a  smoothed  version  of  the  actual  coast  serving 
to  define  one  of  the  coordinate  axes.  The  other  coordinate  can  be 
identified  with  long  wave  travel  time  from  the  seaward  boundary  to  the 
shoreline.  This  representation  provides  good  definition  of  the  open- 
coast  surge,  which  is  the  principal  quantity  sought  in  this  study. 

36.  The  wind  field  model  used  (SPH)  is  not  the  most  theoretically 
advanced  model  currently  available;  but  it  does  account  for  all  impor¬ 
tant  features  of  the  hurricane  wind  field  affecting  the  open-coast  surge, 
and  it  does  so  with  a  very  small  number  of  input  parameters.  Since  this 
particular  study  required  270  different  wind  fields,  economy  of  usage, 
coupled  with  good  accuracy  of  representation,  was  an  important  factor 
in  selection  of  the  SPH  model.  Advanced  planetary  boundary  layer  models 
are  currently  becoming  available  and  they  may  provide  more  accurate  wind 
field  representations,  but  probably  with  increased  complexity  and  cost 
of  implementation.  At  the  present  time,  the  SPH  (with  fetch-dependent 
overland  wind  reduction  factors)  may  be  considered  an  adequate  tool  for 
numerical  stage-frequency  studies. 

37-  After  an  appropriate  wind  field  model  has  been  selected,  one 
is  left  with  the  problem  of  specifying  a  suitable  ensemble  of  wind 
field  parameters.  The  statistics  used  for  definition  of  the  SPH  can  be 
used  with  other  wind  field  models  since  they  are  an  objective  presenta¬ 
tion  of  observed  parameters  within  the  various  coastal  zones  of  the 
Atlantic  Ocean  and  the  Gulf  of  Mexico.  Paucity  of  data  prevents  esti¬ 
mation  of  the  distribution  functions  for  all  pertinent  variables,  lead¬ 
ing  to  specification  of  a  probable  range  as  the  most  reliable  informa¬ 
tion  on  many  of  these  variables.  The  distribution  of  central  pressure 
can  be  estimated,  through  use  of  the  theory  of  extreme-value  statistics, 
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if  a  modest  sized  sample  is  available  (20  or  more  values  may  be 
adequate).  There  is  no  assurance,  however,  that  the  tails  of  the 
distribution  (i.e.,  "unlikely"  events)  can  be  properly  accounted  for. 
The  only  remedy  for  this  problem  is  to  have  a  broader  data  base  from 
which  to  construct  statistical  estimates.  Until  such  data  bases  are 
available,  or  a  substantially  different  methodology  is  introduced, 
the  techniques  used  in  the  present  study  are  believed  to  offer  a  valid 
approach  to  estimation  of  peak  setup  for  open-coast  hurricane  surges. 
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APPENDIX  A:  NOTATION 


Factor  dependent  on  bottom  friction 


BF^.BFgjBF^  Beachfront  sites 

d  Augmented  depth  (mean  depth  plus  projected  peak  surge 
elevation) 

d^  Water  depth  at  breaking  location 

D 

d^  Water  depth  on  seaward  end  of  wave  orthogonal  segment 
di+1  Water  depth  on  landward  end  of  wave  orthogonal  segment 
D  Total  water  depth 
f  Coriolis  parameter 

F  Number  of  peak  surge  or  peak  setup  "Values  less  than  or 
equal  to  W^ 

F,y,v  Variables  scale  factors  describing  the  combined  effects  of 
the  conformal  and  stretching  transformations  applied  to 
the  original  curvilinear  system 

g  Acceleration  due  to  gravity 

H  Sea  surface  elevation  relative  to  a  specified  datum 

IL  Wave  height  on  seaward  end  of  wave  orthogonal  segment 

Hi+1  Modified  wave  height  on  landward  end  of  wave  orthogonal 
segment 

Hg  Hydrostatic  sea  surface  elevation  due  to  atmospheric 
pressure  anomaly  ("inverted  barometer"  effect) 

H,,  Height  of  significant  storm  wave 
b 

H  ,  Combined  peak  surge  and  wave  setup  at  the  beach  (peak 
setup) 


Transfer  function  for  wave  height  across  wave  orthogonal 
segment 


Kg  Linear  long  wave  form  of  the  shoaling  factor 
Wavelength 


A1 


Landfall  sites 


Sample  size  for  estimation  of  exceedance  probabilities 

Number  of  0.5-ft  sampling  windows  required  to  span 
observed  range  of  peak  surge  or  peak  setup 

Cumulative  distribution  function  for  peak  surge  or  peak 
setup 

Volume  transport  per  unit  width 

Radius  to  point  of  maximum  wind  speed 

Peak  surge  at  shoreline 

Setdown  at  the  breaking  zone 

Total,  wave  setup  on  normal  sloping  beach 

Time 

Period  of  significant  storm  wave 
Maximum  wind  speed 

Forward  translational  speed  of  the  storm 

Intrinsic  random  variable  (peak  surge  or  peak  setup) 

Right-hand  value  on  the  nth  sampling  window  for  peak 
surge  or  peak  setup 

Right-hand  value  on  the  Nth  sampling  window 
Length  of  segment  of  wave  orthogonal 

Central  pressure  deficit,  relative  to  constant  peripheral 
pressure 

Wave  setup  between  the  breaking  zone  and  the  shore 
Friction  factor 

Geometric  form  factor  for  shoaling 
Storm  track  heading 

Bottom  stress  (friction)  divided  by  water  density 
Wind  stress  divided  by  water  density 


Shoreline  coordinate;  significant;  shoaling  factor 
Coordinate  normal  to  S  ;  transitional;  track 


In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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